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Alexei Borodin 

Abstract. We introduce and study "isomonodromy" transformations of the matrix 
linear difference equation Y{z + 1) = A{z)Y{z) with polynomial (or rational) A{z). 
Our main result is a construction of an isomonodromy action of — l on the 

space of coefficients A{z) (here m is the size of matrices and n is the degree of A{z)). 
The (birational) action of certain rank n subgroups can be described by difference 
analogs of the classical Schlesinger equations, and we prove that for generic initial 
conditions these difference Schlesinger equations have a unique solution. We also 
show that both the classical Schlesinger equations and the Schlesinger transforma- 
tions known in the isomonodromy theory, can be obtained as limits of our action in 
two different limit regimes. 

Similarly to the continuous case, for m = n = 2 the difference Schlesinger equa- 
tions and their g-analogs yield discrete Painleve equations; examples include dPII, 
dPIV, dPV, and 5-PVI. 



Introduction 

In recent years there has been considerable interest in analyzing a certain class 
of discrete probabilistic models which in appropriate limits converge to well-known 
models of Random Matrix Theory. The sources of these models are quite diverse, 
they include Combinatorics, Representation Theory, Percolation Theory, Random 
Growth Processes, tiling models and others. 

One quantity of interest in both discrete models and their random matrix limits 
is the gap probability - the probability of having no particles in a given set. It is 
known, due to works of many people, see [JMMS], [Me], [TW], [P], [HI], [BD], that 
in the continuous (random matrix type) setup these probabilities can be expressed 
through a solution of an associated isomonodromy problem for a linear system of 
differential equations with rational coefficients. 

The goal of this paper is to develop a general theory of "isomonodromy" trans- 
formations for linear systems of difference equations with rational coefficients. This 
subject is of interest in its own right. As an application of the theory, we show 
in a subsequent publication that the gap probabilities in the discrete models men- 
tioned above are expressible through solutions of isomonodromy problems for such 
systems of difference equations. In the case of one-interval gap probability this has 
been done (in a different language) in [Bor] , [BB] . One example of the probabilistic 
models in question can be found at the end of this Introduction. 

Consider a matrix linear difference equation 

Y{z + l)= A{z)Y{z). (1) 
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Typeset by AmS-^W- 



Here 

A{z) = Aaz'^ + Aiz"-^ + ■ ■ ■ + A^, A, € Mat(m, C), 

is a matrix polynomial and Y : C Mat(m, C) is a matrix meromorphic function.^ 
We assume that the eigenvalues of are nonzero and that their ratios are not real. 
Then, without loss of generality, we may assume that is diagonal. 

It is a fundamental result proved by Birkhoff in 1911, that the equation (1) has 
two canonical meromorphic solutions Y\z) and Y'^{z), which are holomorphic and 
invertiblc for <C and ^ respectively, and whose asymptotics at z = oo 
in any left (right) half-plane has a certain form. Birkhoff further showed that the 
ratio 

P{z) = {Y^{z))-'Y\z), 

which must be periodic for obvious reasons, is, in fact, a rational function in 
exp(27rz2;). This rational function has just as many constants involved as there 
are matrix elements in ^i, . . . , j4„. Let us call P{z) the monodromy matrix of (1). 

Other results of Birkhoff show that for any periodic matrix P of a specific form, 
there exists an equation of the form (1) with prescribed Aq, which has P as the 
monodromy matrix. Furthermore, if two equations with coefficients A{z) and A{z), 
Aq = Aq, have the same monodromy matrix, then there exists a rational matrix 
R{z) such that 

Aiz) = Riz + l)Aiz)R~\z). (2) 

The first result of this paper is a construction, for generic A{z), of a homomor- 
phism of X'raC^+i)-! into the group of invcrtible rational matrix functions, such that 
the transformation (2) for any R{z) in the image, does not change the monodromy 
matrix. 

If we denote by ai , . . . , Umn the roots of the equation det A{z) =0 (called eigen- 
values of A{z)) and by di, . . . , d„ certain uniquely defined exponents of the asymp- 
totic behavior of a canonical solution Y{z) of (1) at z = oo, then the action of 
jm{n+i)-i jg uniquely defined by integral shifts of {a^} and {dj} with the total 
sum of all shifts equal to zero. (We assume that Ui — aj ^ Z and di — dj ^ Z for 
any i 7^ j.) 

The matrices R{z) depend rationally on the matrix elements of {Ai}"^]^ and 
{oj}^" (^0 is always invariant), and define birational transformations of the vari- 
eties of {Aj} with given {ui} and {dj}. 

There exist remarkable subgroups Z" C z^^^+i)-! which define birational trans- 
formations on the space of all A{z) (with fixed Aq and with no restrictions on the 
roots of det A{z)), but to see that we need to parameterize A{z) differently. 

To define the new coordinates, we split the eigenvalues of A{z) into n groups of 
m numbers each: 

{ai,...,a„„} = {a(^\...,a(^)}U---U{4"\...,aW}. 
The splitting may be arbitrary. Then we define Bi to be the uniquely determined 

such that z — Bi is a right divisor of A{z): 

A{z) = {Aoz^-^ + A[z"-^ + ■■■ + A'^_,){z - Bi). 



^Changing Y{z) to (T{z))''Y(z) readily reduces a rational A.(z) to a polynomial one. 
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The matrix elements of {Bi}f^^ are the new coordinates on the space of A(z). 

The action of the subgroup Z" mentioned above consists of shifting the eigen- 
values in any group by the same integer assigned to this group, and also shifting 
the exponents {d;} by the same integer (which is equal to minus the sum of the 
group shifts). If we denote by {Bi{ki, . . . , kn)} the result of applying fc € Z" to 
{Bi}, then the following equations are satisfied: 

Bi{...)-B,{...,k, + l,...) = B,i...)-B,i...,h + l,...), (3) 
Bj{. ..,ki + l,... )B,{. . . ) = B,{. ..,% + !,... )B,{. . . ), (4) 
Bi{ki + 1, . . . , fc„ + 1) = A^^Bi{ki, kn)Ao - I, (5) 

where i,j = l,...,n, and dots in the arguments mean that other fcj's remain 
unchanged. We call them the difference Schlesinger equation for the reasons that 
will be clarified below. Note that (3) and (4) can be rewritten as 

{z-Bi{...,kj + l,...)){z-Bj{...)) = {z-Bj{...,ki + l,...)){z-Bi{...)). 

Independently of Birkhoff 's general theory, we prove that the difference Schlesinger 
equations have a unique solution satisfying 

Sp{Bi{ki,...,kn)) = Sp{Bi)-ki, i = l,...,n, (6) 

for an arbitrary nondegenerate Aq and generic initial conditions {Bi — Bi{0)}. 
(The notation means that the eigenvalues of Bi{k) are equal to those of Bi shifted 
by —/;•;.) Moreover, the matrix elements of this solution are rational functions in 
the matrix elements of the initial conditions. This is our second result. 

In order to prove this claim, we introduce yet another set of coordinates on A{z) 
with fixed ^o. which is related to {Bi} by a birational transformation. It consists 
of matrices Ci G Mat(m,C) with Sp{Ci) = Sp{Bi) such that 

A{z) = Ao{z-C^)---{z-Cn). 

In these coordinates, the action of Z" is described by the relations 

{z + l-Ci)---{z + l- C,,)Aa{z - Ci) • • • (z - C,_i) 
= {z + l- Q+i) • • • (z + 1 - C„)Ao(z -Ci)---{z- a), (7) 
Cj = Cj{ki, kn), Cj = Cj{ki, ki-i,ki + 1, ki+i, ...,kn) for all j. 

Again, we prove that there exists a unique solution to these equations satisfying 
Sp{Ci{k)) = Sp{Ci) — ki, for an arbitrary invertible Aq and generic {Ci = Ci{0)}. 
The solution is rational in the matrix elements of the initial conditions. 

The difference Schlesinger equations have an autonomous limit which consists of 

(3), (4), and 

Bi{ki + 1,.. .,fc„ + 1) = AQ^Bi{ki,...,kn)Ao, (5-aut) 

Sp{Bi{ki,...,kn)) = Sp{Bi), i = l,...,n. (6-aut) 
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The equation (7) then turns into 



{z-Ci)---{z^ Cn)Ao{z - Ci) • • • (z - 
{z - Q+i) • • • (z - Cn)Ao{z - Ci) • • • (z - C\). 



(7-aut) 



The solutions of these equations were essentially obtained in [V] via a general 
construction of commuting flows associated with set-theoretical solutions of the 
quantum Yang-Baxter equation, see [V] for details and references. 

The autonomous equations can also be explicitly solved in terms of abelian func- 
tions associated with the spectral curve {{z,w) : det{A{z) — wl) — 0},^ very much 
in the spirit of [MV, §1.5]. We hope to explain the details in a separate publication. 

The whole subject bears a strong similarity (and not just by name!) to the 
theory of isomonodromy deformations of linear systems of differential equations 
with rational coefBcients: 



which was developed by Schlesinger around 1912 and generalized by Jimbo, Miwa, 
and Ueno in [JMU], [JM] to the case of higher order singularities. If we analytically 
continue any fixed (say, normalized at a given point) solution 3^(C) of (8) along a 
closed path 7 in C avoiding the singular points {xk} then the columns of y will 
change into their linear combinations: y 1— > yM.y. Here is a constant invertible 
matrix which depends only on the homotopy class of 7. It is called the monodrom,y 
matrix corresponding to 7. The monodromy matrices define a linear representation 
of the fundamental group of C with n punctures. The basic isomonodromy problem 
is to change the differential equation (8) so that the monodromy representation 
remains invariant. 

There exist isomonodromy deformations of two types: continuous ones, when 

Xi move in the complex plane and Bi = Bi{x) form a solution of a system of 
partial differential equations called Schlesinger equations, and discrete ones (called 
Schlesinger transformations), which shift the eigenvalues of Bi and exponents of 
3^(^) at C = CO by integers with the total sum of shifts equal to 0. 
We prove that in the limit when 



our action of 2'"'"+^^"^ in the discrete case converges to the action of Schlesinger 
transformations on Bj. This is our third result. 

Furthermore, we argue that the "long-time" asymptotics of the Z"-action in 
the discrete case (that is, the asymptotics of Bi{[xie~^], . . . , [a;„e~-^])), e small, is 
described by the corresponding solution of the Schlesinger equations. More exactly, 
we conjecture that the following is true. 

Take Bi = Bi{e) G Mat(m, C), i = 1, . . . ,n, such that 





Bi =Xie ^ + Bi 



'11 



e ^ 0, 



Bi{e) - ViC 



-1 



+ Bi^Q, 



e 



0. 



^It is easy to see that the curve is invariant under the flows. 

4 



Let Bi{ki, . . . , fc„) bo the solution of the difference Schlesinger equations (3.1)-(3.3) 
with the initial conditions {-Bi(O) = Bi}, and let Bi{xi, . . . be the solution of 
the classical Schlesinger equations (5.4) with the initial conditions . . . , 2/„) = 

Bi}. Then for any xi,. . . ,Xn € K and i = 1, . . . , n, we have 

Bi([xie~'^], . . . , [x„e~'^]) + (xi - yi)e~^ + Bi{yi -xi,.. . , j/„ - a;„) ^ 0, e ^ 0. 

In support of this conjecture, we explicitly show that the difference Schlesinger 
equations converge to the conventional Schlesinger equations in the limit e ^ 0. 

Note that the monodroniy representation of 7ri(C \ {xi, . . . , Xn}) which provides 
the integrals of motion for the Schlesinger flows, has no obvious analog in the dis- 
crete situation. On the other hand, the obvious differential analog of the periodic 
matrix P, wliic;li contains all integrals of motion in the case of difference equa- 
tions, gives only the monodromy information at infinity and does not carry any 
information about local monodromies around the poles Xi, . . . ,Xn- 

Most of the results of the present paper can be carried over to the case of q- 
difference equations of the form Y{qz) = A{z)Y{z). The g-difference Schlesinger 
equations are, cf. (3-6), 

Bi{...)-Bi{...,q>'^+\...)=Bj{...)-Bj{...,q''^+\...), (3q) 

Bj{. . . , q^^+\ . . . )B,{. ..) = Bi{. . . , g'^^+i, . . . )Bj{. . . ), (4q) 

B.(g'=^+\ . . . , g'^^+i) = q-'A^'B,{q'^\. . . , g'=")Ao, (5q) 

Sp{Bi{q''\...,q^-))=q-'''Sp{Bi), i = l,...,n. (6q) 

The g-analog of (7) takes the form 

{z - q-^Ci) • • • (z - g-iC„)Ao(z - Ci) • • • (z - Q^i) 
= {z-q-^C,+i)---{z-q-^Cr)Ao{z-Ci)---{z~C,), (7q) 
Cj = Cj {q''^ ,...,q'-), Cj = Cj {q"^ q"-^ , q'''+' , , . . . , ) for all j. 

A more detailed exposition of the ^'-difference case will appear elsewhere. 

Similarly to the classical case, see [JM], discrete Painleve equations of [JS], 
[Sak] can be obtained as reductions of the difference and g-difference Schlesinger 
equations when both m (the size of matrices) and n (the degree of the polynomial 
A{z)) are equal to two. For examples of such reductions see [Bor, §3] for difference 
Painleve II equation (dPII), [Bor, §6] and [BB, §9] for dPIV and dPV, and [BB, 
§10] for q-PYl. This subject still remains to be thoroughly studied. 

As was mentioned before, the difference and g-difference Schlesinger equations 
can be used to compute the gap probabilities for certain probabilistic models. We 
conclude this Introduction by giving an example of such a model. We define the 
Hahn orthogonal polynomial ensemble as a probability measure on all Z-point subsets 
of {0,1,...,N}, N > I > 0, such that 

I 

Prob{(a;i, ... ,a;/)} = const- JJ {xi - xj)"^ ■ Y\_u!{x^), 

l<i<j<l i=l 
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where w{x) is the weight function for the classical Hahn orthogonal polynomials: 
w{x)={ ,r ' q;,/3 > -1 or a,/3 < -A/^. 



X J \ N — X 

This ensemble came up recently in harmonic analysis on the infinite-dimensional 
unitary group [BO, §11] and in a statistical description of tilings of a hexagon by 
rhombi [Joh, §4]. 

The quantity of interest is the probability that the point configuration {xi, . . . ,Xi) 
does not intersect a disjoint union of intervals [/ci, A:2] LJ • • • U [fc2s-i, fc2s]- As a func- 
tion in the endpoints k\,. . . ,k2s G {0, 1, ... , A^}, this probability can be expressed 
through a solution of the difference Schlesingcr equations (3)- (6) for 2 x 2 matrices 
with n = degA{z) = s + 2, Aq = /, 

Sp{B.i) = {-ki,-ki}, i = l,...,2s, 
Sp{B2s+i) U Sp{B2s+2) = {Q,-a,N+l,N + l + /?}, 

and with certain explicit initial conditions. The equations are also suitable for 
numerical computations, and we refer to [BB, §12] for examples of those in the case 
of one interval gap. 

I am very grateful to P. Deift, P. Deligne, B. Dubrovin, A. Its, D. Kazhdan, 
I. Krichever, G. Olshanski, V. Retakh, and A. Veselov for interesting and helpful 
discussions. 

This research was partially conducted during the period the author served as a 
Clay Mathematics Institute Long-Term Prize Fellow. 

1. Birkhoff's theory 

Consider a matrix linear difference equation of the first order 

Y{z + l) = A{z)Y{z). (1.1) 

Here A : <C ^ Mat(m, C) is a rational function (i.e., all matrix elements of A{z) 
are rational functions of z) and m > 1. We are interested in matrix meromorphic 
solutions y : C — > Mat(m, C) of this equation. 

Let n be the order of the pole of A{z) at infinity, that is, 

A{z) = Aqz'" + Aiz"~^ + lower order terms . 

We assume that (1.1) has a formal solution of the form 

Y{z) = ^"^e-"^ (yo + ^ + ^ + ...^ diag [p^ z''\ . . . , p^, z''-) (1.2) 

with pi,. . . ,pm ^0 and det Yq ^ 0.^ 

It is easy to see that if such a formal solution exists then pi,. . . ,pm must be the 
eigenvalues of Aq, and the columns of % must be the corresponding eigenvectors 
of Aq. 

Note that for any invertible T G Mat(m, C), {TY){z) solves the equation 

{TY){z + 1) = {TA{z)T-^) {TY){z). 

Thus, if Aq is diagonalizable, we may assume that it is diagonal without loss of 
generality. Similarly, if Aq = I and Ai is diagonalizable, we may assume that Ai 
is diagonal. 

^Substituting (1.2) in (1.1) we use the expansion (^)"^ = e"^'°(i+^"') = e" - ^ + . . . 
to compare the two sides. 
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Proposition 1.1. If Aq = diag(/7i, . . . ,Pm), where {pi}^i are nonzero and pair- 
wise distinct, then there exists a unique formal solution of (1.1) of the form (1.2) 
mth Yq = I. 

Proof. It suffice to consider the case n = 0; the general case is reduced to it by 
considering {r{z))"'Y{z) instead ofY{z), because 

r(^) = \/2^^^-ie-^ {^ + ^^~^ + ---^- 

(More precisely, this expression formally solves T{z + 1) = zT{z).) 

Thus, we assume n = 0. Then we substitute (1.2) into (1.1) and compute Yk one 
by one by equating the coefficients of z^', Z = 0, 1, . . . . If lb = / then the constant 
coefficients of both sides are trivially equal. The coefficients of z~^ give 

YiAq + diag(pirfi, . . .,Pmdm) = ^o^i + ^i- (1-3) 

This equality uniquely determines {di} and the off-diagonal entries of Yi, because 

[Yi,Ao]^j = {pj-pi){Yi),j. 

Comparing the coefficients of z~'^ we obtain 

(F2 - Yi)Ao + Yi diag(pidi, . . . , pmdm) + • • • = A0Y2 + A^Yi + ... , 

where the dots stand for the terms which we already know (that is, those which 
depend only on p^'s, rfj's, A^'s, and Yq ^ I). Since the diagonal values of Ai are 
exactly pidi, . . . Pndn by (1.3), we see that we can uniquely determine the diagonal 
elements of li and the off-diagonal elements of Y2 from the last equality. 

Now let us assume that wc already determined Yi, . . . , J/_2 and the off-diagonal 
entries of by satisfying (1.1) up to order I — 1. Then comparing the coefficients 
of z~^ we obtain 

{Yi -{I- 1)1^-1)^0 + Yi-i diag(/9idi, . . . , p^d^) + . . . = A^Yi + AiY-i + • • • , 

where the dots denote the terms depending only on pj's, dj's, Aj's, and lb) • • • > li-2- 
This equality allows us to compute the diagonal entries of YJ-i and the off-diagonal 
entries of Yi . Induction on I completes the proof. □ 

The condition that the eigenvalues of Aq arc distinct is not necessary for the 
existence of the asymptotic solution, as our next proposition shows. 

Proposition 1.2. Assume that Aq = I and A\ — diag(ri, . . . , r„) where ri — rj ^ 
{±1, ±2, . . . } for all i,j = l,...,n. Then there exists a unique formal solution of 
(1.1) of the form (1.2) with Yq = I. 

Proof. As in the proof of Proposition 1.1, we may assume that n = 0. Comparing 
constant coefficients we see that pi = • • • = = 1- Then equating the coefficients 
of z~^ we find that di = ri, i = 1, . . . ,m. Furthermore, equating the coefficients of 
2;"', Z > 2 we find that 

[Yi-uAi]-{l-l)Yi-i 

is expressible in terms of Ai's and Yi,. . . , Yi-2. This allows us to compute all l^'s 
recursively. □ 

We call two complex numbers zi and Z2 congruent if zi — Z2 & Z. 

7 



Theorem 1.3 (G. D. BirkhofF [Bil, Theorem III]). Assume that 



Ao = diag(pi, . . .,Prn), 

Pi ^ 0, i = 1, . . . , TO, Pi/pj for all i =^ j. 

Then there exist unique solutions Y^z) {Y'^{z)) of (1.1) such that 

(a) the function y'(z) {Y^{z)) is analytic throughout the complex plane except 
possibly for poles to the right (left) of and congruent to the poles of A{z) {respec- 
tively, A^^{z - 1)); 

(b) in any left {right) half plane Y''{z) {Y^{z)) is asymptotically represented by 
the right-hand side of (1.2). 

Remeirk 1.4. Part (b) of the theorem means that for any fc = 0, 1, . . . , 



Y^'^{z) ^-"^e"^ diag(pr^^"'^' , ■ 



vfe-l 



< 



const 



yk 



for large in the corresponding domain. 

Theorem holds for any (fixed) choices of branches of ln(z) in the left and right 
half planes for evaluating z~"^ = e~"^'"(^) and z~'^'' — e~'''='"(^\ and of a branch 
of ln(p) with a cut not passing through pi,. . . ,pm for evaluating p^^ = e~^''^'"=. 
Changing these branches yields the multiplication o{Y^''^{z) by a diagonal periodic 
matrix on the right. 

Remark 1.5. Birkhoff states Theorem 1.3 under a more general assumption: he 
only assumes that the equation (1.1) has a formal solution of the form (1.2). How- 
ever, as pointed out by P. Deligne, Birkhoff 's proof has a flaw in case one of the 
ratios pi / pj is real. The following counterexample was kindly communicated to me 
by Prof. Deligne. 

Consider the equation (1.1) with m = 2 and 



A{z) = 

The formal solution (1.2) has the form 



1 1/z 
1/e 



z-' + 



1 

e-' 



with a — e/(l — e). 

Actual solutions that we care about have the form 



Y{z) = 



1 u{z) 
e-^ 



where u{z) is a solution of u{z + 1) = u{z) + e ^/z. In a right half-plane we can 
take 



u^{z) 



ra=0 



Z + n 



e ae 



The first order approximation of u^{z) anywhere except near nonpositive integers 
is 

°° p-(z+n) 

Next terms can be obtained by expanding l/{z + n). 

In order to obtain a solution which behaves well on the left, it suffices to cancel 
the poles: 

u\z)=u'^{z) 

The corresponding solution Y^{z) has needed asymptotics in sectors of the form 
7r/2 + e < argz < 37r/2 + e, but it has wrong asymptotic behavior as z ^ +ioo. 
Indeed, lim^^+joo u^z) = —2m. 
On the other hand, we can take 

u\z) = u\z) + 2m = u'{z) + ^——^ , 

which has the correct asymptotic behavior in 7r/2 — e < axgz < 3^/2 — e, but fails 
to have the needed asymptotics at —ioo. 

Remark 1.6. In the case when \pi\ > \p2\ > ■ ■ ■ > \pm\ > 0, a result similar to 
Theorem 1.3 was independently proved by R. D. Carmichael [C]. He considered 
the asymptotics of solutions along lines parallel to the real axis only. Birkhoff also 
referred to [N] and [G] where similar results had been proved somewhat earlier. 

Now let us restrict ourselves to the case when A{z) is a polynomial in z. The 
general case of rational A{z) is reduced to the polynomial case by the following 
transformation. If {z — Xi) ■ ■ ■ {z — Xg) is the common denominator of {Ahi{z)} (the 
matrix elements of A{z)), then 

Y{z) = T{z - xi) ■ ■ - Viz - Xs) -Yiz) 

solves Y{z + 1) = A{z)Y{z) with polynomial 

A{z) = {z-xi)---{z- Xs)A{z). 

Note that the ratio P{z) = {Y^{z))~^ Y\z) is a periodic function. (The relation 
P{z + 1) = P(z) immediately follows from the fact that F'^'' solve (1.1).) From 
now on let us fix the branches of ln(2;) in the left and right half planes mentioned 
in Remark 1.4 so that they coincide in the upper half plane. Then the structure of 
P{z) can be described more precisely. 

Theorem 1.7 (G. D. Birkhoff [Bil, Theorem IV]). In the assumptions of The- 
orem 1.3, the matrix elements pki{z) of the periodic matrix P{z) = {y^ {z))~^Y'' {z) 
have the form 

Pkkiz) = 1 + 4^e2-- + . . . + c^-^-i^e^-^"-!)- + e2-'^'=e2.™.^ 
p^i{z) = 6^-^'='^ (4°) + 4f e2-« + . . . + c(7-i)e2-(-i)-) {k + I), 



where c^^ are some constants, and Xki denotes the least integer as great as the real 
part of (ln(/9/) - \n{pk))/2ni. 

Thus, starting with a matrix polynomial A{z) — Aqz^ + Aiz"'~^ + • • • + A„ 
with nondegenerate Aq = diag(/Oi, . . . , pm), pk ^ pi for k I, we construct the 
characteristic constants {dk}, {^4^'} using Proposition 1.1 and Theorems 1.3, 1.7. 

Note that the total number of characteristic constants is exactly the same as 
the number of matrix elements in matrices Ai, . . . , An- Thus, it is natural to ask 
whether the map 

is injective or surjectivc (the constants pi, . . . ,/>„ are being fixed). The following 
partial results are available. 

Theorem 1.8 (G. D. Birkhoff [Bi2, §17]). For any nonzero pi,. ■ . ,pm, Pi/Pj ^ 
M for i ^ j, there exist matrices Ai, . . . ,An such that the equation (1.1) with Aq = 
diag{pi, ... pm) either possesses prescribed characteristic constants {dk} , {c\^i}, or 
else constants {dk + h}, {c\a'}> where li,...,lrn o-re integers. 

Theorem 1.9 (G. D. Birkhoff [Bil, Theorem VII]). Assume that we have two 
matrix polynomials A\z) = A'qz"^ -| \- A'^ and A"{z) = A!^z'^ ^ h with 

^0 = K = diag(pi, . . . Pm), Ph 7^ 0, pujpi ^^fork^l, 

such that the sets of the characteristic constants for the equations Y'{z + 1) = 
A'{z)Y'{z) and Y"{z + 1) = A" {z)Y" {z) are equal. Then there exists a rational 
matrix R{z) such that 

A"{z) = R{z + l)A'{z)R-\z), 

and the left and right canonical solutions Y^''^ of the second equation can he obtained 
from those of the first equation by multiplication by R on the left: 

(Y"\hr ^ D (Y'\l,r 



2. ISOMONODROMY TRANSFORMATIONS 

The goal of this section is to construct explicitly, for given A{z), rational matrices 
R{z) such that the transformation A{z) R{z + \)A(z)R^^(z), cf. Theorem 1.9 
above, preserves the characteristic constants (more generally, preserves {c^'^'} and 
shifts dfc's by integers). 

Let A{z) be a matrix polynomial of degree n> 1^ A^ = diag(pi, . . . , Pm), and pi's 
are nonzero and their ratios are not real. Fix mn complex numbers ai , . . . , Omn 
such that ai — aj ^ % for any i ^ j. Denote by M.{a\,. . . , Om„; di, . . . , rf^) the 
algebraic variety of all n-tuples of m by m matrices Ai, . . . , A„ such that the scalar 
polynomial 

det A{z) = det(ylo-2" + ^i-z""^ + • • • + ^„) 

of degree mn has roots ai, . . . , amn-, and pi [di — f ) = {.A\)^^ (this comes from the 
analog of (1.3) for arbitrary n). 

10 



Theorem 2.1. For any /^i, . . . , Kmn € ^i, . . . , ^ ^? 

mn 7n 

there exists a nonempty Zariski open subset A of M.{ai, . . . , Omn', di, . . . ,dm) such 
that for any {Ai, . . . , A„) e A there exists a unique rational matrix R{z) with the 
following properties: 

A{z) = R{z + l)A{z)R-\z) = lo^" + ^1-2""' + ••• + !„, Ao = Ao, 

{Ai, . . .,An) e M{ai +Ki,.. .,amn + Kmn;dl +Si,...,dm + Sm), 

and the left and right canonical solutions of Y{z + 1) = A{z)Y{z) have the form 

where Y^'"'' are left and right canonical solutions ofY{z + 1) = A{z)Y{z). 

The map {Ai, . . . , An) f-^ (Ai, . . . , An) is a birational map of algebraic varieties. 

Remark 2.2. The theorem implies that the characteristic constants {c^'^'} for the 
difference equations with coefficients A and A are the same, while the constants dk 
are being shifted by 6k € 

Note also that if we require that all dk 's do not change then, by virtue of The- 
orem 1.9, Theorem 2.1 provides all possible transformations which preserve the 
characteristic constants. Indeed, if A"{z) = R{z + 1)A' {z)R~^{z) then zeros of 
det^"(z) must be equal to those of dety4'(z) shifted by integers. 

Proof. Let us proof the uniqueness of R first. Assume that there exist two rational 
matrices Ri and R'2 with needed properties. This means, in particular, that the 
determinants of the matrices 

1(1) =Ri{z + l)A{z)R^\z) and I^^) = r^[z 

vanish at the same set of mn points Oi = Oi + Ki, none of which are different by 
an integer. Denote by Y{ = R^Y^ and Y2 = R2Y^ the right canonical solutions of 
the corresponding equations. Then Y{{Y2)~^ = iJii?^""^ is a rational matrix which 
tends to / at infinity. Moreover, 

{Y{i%l-') {z + l)= A('\z) {y{{Y2^)-') {z) • 

Since Y{(Y{)^^ is holomorphic for 3?z ^ 0, the equation above implies that this 
function may only have poles at the points which are congruent to Oi (zeros of 
detA^'^^z)) and to the right of them. (Recall that two complex numbers are con- 
gruent if their difference is an integer.) But since (YJ)"^ is also holomorphic for 
?Siz ^ 0, the same equation rewritten as 

{y{in)-') (z) = (aW(^))"' {y{{%^)-') {z + 1)A('\z) 

implies that this fimction may only have poles at the points (zeros of dct A^-'^^ (2;)) 
or at the points congruent to them and to the left of them. Thus, Y{{Y{)~^ = 
RiR^^ is entire, and by Liouville's theorem it is identically equal to /. The proof 
of uniqueness is complete. 

To prove the existence we note, first of all, that it suffices to provide a proof if 
one of Kj's is equal to ±1 and one of Sj's is equal to =f1 with all other k's and ^'s 
being equal to zero. The proof will consist of several steps. 
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Lemma 2.3. Let A{z) he an m by m matrix-valued function holomorphic near 
z = a, and detA{z) = c{z — a) + O ({z — a)^) as z a, where c 0. Then there 
exists a unique (up to a constant) nonzero vector v G C" such that A{a)v = 0. 
Furthermore, if B{z) is another matrix valued function which is holomorphic near 
z = a, then {BA^^){z) is holomorphic near z = a if and only if B(a)v = 0. 

Proof. Let us denote by Ei the matrix unit which has 1 as its (1, l)-entry and as 
all other entries. Since det^(a) = 0, there exists a nondegenerate constant matrix 
C such that A{a)CEi = (the first column of C must be a 0-eigenvector of A{a)). 
This implies that 

H{z) = A{z)C{Ei{z - a)-^ + 1 - Ei) 

is holomorphic near z = a. On the other hand, detH{a) = cdetC ^ 0. Thus, 
A(a) = H(a){I — Ei)C^^ annihilate a vector v if and only if C^^v is proportional 
to (1,0,..., 0)*. Hence, v must be proportional to the first column of C. The proof 
of the first part of the lemma is complete. 
To prove the second part, we notice that 

{BA-^){z) = B{z)C{Ei{z - a)-i + 7 - Ei)H-\z) 

which is bounded at z = a if and only if B{a)CEi =0. □ 

More generally, we will denote by Ei the matrix unit defined by 



{Ei) 



kl 



1, k = I ^ i, 
0, otherwise. 



Lemma 2.4 ([JM, §2 and Appendix A]). For any nonzero vector v — (wi, . . . , VmY , 
Q € Mat(TO, C), a € C, and i € {l,...,m}, there exists a linear matrix-valued 
function R{z) = R-i{z — a) + Rq with the properties 

R{z) (I + Qz-^ + O (0-2)) z-^^ =1 + 0{z-^), ^ oo, 
Rov = 0, 

if and only ifvi^O. In this case, R'^^{z) is given by 

R{z) = E,{z - a) + Ro, R^^ {z) = I - E, + Ri{z - aY^ , deiR{z) = z-a, 

"^s^i Qis'^^st k = I = i, 

~Qii^ k = I if 

-v~^Vk, k^i, l^i, 

. Skh k ^i, I ^ i, 

R\ = V^ [Qili • ■ • ) Qi,i-1, 1, Qi,i+1, ■ ■ ■ ■> Qim ] ■ 



(Ro) 



kl 



The proof is straightforward. 

Now wc return to the proof of Theorem 2.1. 

Assume that ki = — 1, 5i = 1 for some i = 1, . . . , m, and all other k's and J's are 
zero. Since ai is a simple root of det^(2;), by Lemma 2.3 there exists a unique (up 
to a constant) vector v such that A{a)v = 0. Clearly, the condition Vi ^ defines 
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a nonempty Zariski open subset of A4{ai, . . . , Qmn', <5i, • • • ; Sm)- On this subset, let 
us take R{z) to be the matrix afforded by Lemma 2.4 with o = Oi and Q = Yi 
(we assume that Yq = I, see Proposition 1.1). Then by the second part of Lemma 
2.3, {A{z)R~^{z))~^ = R{z)A~^{z) is holomorphic and invertible near z = ai (the 
invertibility follows from the fact that dot i?(z)^~^(z) tends to a nonzero value as 
z ai). Thus, A{z) = R{z + l)A{z)R^^{z) is entire, hence, it is a polynomial. 
Since 

det A(2;) = — — detA{z) = c{z + 1 ~ ai){z — 02) ■ ■ ■ {z — amn), 0, 
z — ai 

the degree of A{z) is > n. Looking at the asymptotics at infinity, we see that 
deg A{z) < n, which means that A is a polynomial of degree n: 

A{z) = Aoz" + ••• + !„, Ao^O. 

Denote by y'^'' the left and right canonical solutions of y(z + l) = A{z)Y{z) (see 
Theorem 1.3 above). Then := RY^''' are solutions of Y{z + 1) = A{z)Y{z). 
Moreover, their asymptotics at infinity at any left (right) half plane, by Lemma 2.4, 

is given by an expansion of the form (1.2) with Yq = I , p^. = pi^ for all A; = 1, . . . , m, 
and 

~ ( dk + 1, k = i, 
I dk, kj^ i. 

This implies that Aq = diag(pi, . . . , pm), and that F'''" are the left and right canon- 
ical solutions of the equation Y{z + 1) = A{z)Y{z). Indeed, their asymptotic 
expansion at infinity must also be a formal solution of the equation, the fact that 
are holomorphic for JR^; <C {'^ 0) follows from the analogous property for 
F'"'", and the location of possible poles of y is easily determined from the equation. 

For a future reference let us also find a (unique up to a constant) vector v such 
that A*(ai-l)t; = 0. This means that i?~*(ai - l)A*(ai - l)i?*(ai) z; = 0. Lemma 
2.4 then implies that 

V = [(Fi),i, . . . , (Fi). 1, (Fi). ,^1, . . . , (Fi),^]* 

is a solution. Note that Vi 7^ 0. 

Now let us assume that ki = 1 and 5i = —1 for some i = 1, . . . ,m. By Lemma 
2.3, there exists a unique (up a to a constant) vector w such that A*'{ai)w = 0. The 
condition Wi ^ defines a nonempty Zariski open subset of A^(ai, . . . , amn', ^1, • ■ • Sm) 
On this subset, denote by R'{z) the rational matrix-valued function afforded by 
Lemma 2.4 with a = ai, v = w, and Q = — F/ (again, we assume that Fq = /). Set 

R{z) := {Rr\z-1). 

Then by Lemma 2.4 

R{z) (j + Yiz-^ + O (^-2)) z^' =1 + 0{z-^), z^oo. 
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Furthermore, by Lemma 2.3, R *(z + 1)A *{z) is holomorphic and invcrtiblc near 
z = ai. Hence, A{z) = R{z + l)A{z)R~^{z) is entire (note that R~^{z) = {R'f{z- 
1) is Hnear in z). The rest of the argument is similar to the case k\ = —1, 5i = 1 
considered above. ^ 

Finding a solution w to A{a\ + l)u; = is equivalent to finding a solution to 
R'{a\)w = 0. One such solution has the form 

^ = [-iyi)u^ ■ ■ ■ , 1, • • • ' -(^^i)™ ]* 

and all others are proportional to it. Note that its iih coordinate is nonzero. 
From what was said above, it is obvious that the image of the map 

M{ai, .. ., amn] Si,. . . , Sm) M{ai - 1, . . .,amn; Si, . . . ,Si + 1, . . .,Sm) 

is in the domain of definition of the map 

M{ai - 1, . . . , amn, Si,...,S^ + l,..., S^) M{ai, a„„; Si,...,Sm) 

and the other way around. On the other hand, the composition of these maps in 
either order must be equal to the identity map due to the uniqueness argument 
in the beginning of the proof. Hence, these maps are inverse to each other, and 
they establish a bijection between their domains of definition. The rationality of 
the maps follows from the explicit formula for R{z) in Lemma 2.4. The proof of 
Theorem 2.1 is complete. □ 

Remark 2.5. Quite similarly to Lemma 2.4, the multiplier R{z) can be computed 
in the cases when two k's are equal to ±1 or two S's are equal to ±1 with all other 
k's and ^'s being zero, cf. [JM]. 

Assume = — 1 and Kj ~ 1. Denote by v and w the solutions of A{ai) v = 
and A*{aj) w = 0. Then R exists if and only if {v, w) := v*w = w^v ^ 0, in which 
case 

R{z) = I + R-\z) = I-^^, det R{z)= 

z — Gj — 1 z — ai z — aj — 1 

Ro = —, ^ — vvf. 

{v,w) 

Now assume Si = 1, Sj = —1. Then wc must have det R{z) = 1 and 

R{z) (l + Yiz-^ + Y2Z-^ + 0(^-3)) z^'-^' =1 + 0{z-^), z^oo. 

The solution exists if and only if (Yi)-- ^ 0, in which case it has the form 

R{z) = E,z + Ro, R-\z) = EjZ + Ro\ 
with {Ro)ki given by 

l = i l^j 

k^ij Ski 
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and {Rq ^)^( given by 



k 



1 = 1 









k i,j 



k 



J 







1 



(ilk 



3. Difference Schlesinger equations 



In this section we give a different description for the transformations A i— > A of 
Theorem 2.1 with 



and all other equal to zero, and for compositions of such transformations. 

In what follows wc always assume that our matrix polynomials A{z) = Aqz"' + . . . 
have nondegenerate highest coefficients: det Aq ^ 0. We also assume that mn roots 
of the equation det A{z) = are pairwise distinct; we will call them the eigenvalues 
of A(z). For an eigenvalue a, there exists a (unique) nonzero vector v such that 
A{a) V = 0, see Lemma 2.3. We will call v the eigenvector of A{z) corresponding 
to the eigenvalue a. The word generic everywhere below stands for "belonging to 
a Zariski open subset" of the corresponding algebraic variety. 

We start with few simple preliminary lemmas. 

Lemma 3.1. The sets of eigenvalues and corresponding eigenvectors define A[z) 
up to a multiplication by a constant nondegenerate matrix on the left. 

Proof. If there are two matrix polynomials A' and A" with same eigenvalues and 

eigenvectors, then {A' {z))~^ A" (z) has no singularities in the finite plane. Moreover, 
since the degrees of A'{z) and A"{z) are equal, {A'{z))~'^A"{z) ~ {A'q)-^A'^ as 
.2 ^ DO. Liouville's theorem concludes the proof. □ 

We will say that z — B, B G Mat(m, C), is a right divisor of A{z) if A{z) = 
A{z){z — B), where A{z) is a polynomial of degree n — 1. 

Lemma 3.2. A linear function z — B is a right divisor of A{z) if and only if 



Proof See, e.g., [GLR]. □ 

Lemma 3.3. Let ai, . . . , am be eigenvalues of A{z) and vi, . . . , Vm be the cor- 
responding eigenvectors. Assume that vi, . . . ,Vm are linearly independent. Take 
B e Mat(m, C) such that Bvi = aiVi, i = 1, . . . ,m. Then z — B is a right divisor 
of A{z). Moreover, B is uniquely defined by the conditions that z — B is a right 
divisor of A{z) and Sp{B) = {ai, . . . ,am}- 

Proof. For alH = 1, . . . , m, we have 



{AoB" + A^B^-^ + ■■■+ A^)vi = (Ao< + + • • • + ^n)vi = A{ai)vi = 0. 
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Kii = • • • = Kin, = h = ■ ■ ■ = = Tl, 



AqB^ + + • • • + A„ = 0. 



Lemma 3.2 shows that z — B is a right divisor of A{z). 
To show uniqueness, assume that 

A{z) = A'{z){z - B') = A"{z){z - B"). 

This imphes {A" {z))-'^ A' {z) = [z - B"){z - B')"^- Possible singularities of the 
right hand side of this equality are z = ai, i = 1, .... m, while possible singularities 
of the left hand side are all other eigenvalues of A{z). Since the eigenvalues of A(z) 
are pairwise distinet, both sides are entire. But (z — B"){z — B')~^ tends to I as 
z 00. Hence, by Liouville's theorem, B' = B" . □ 

Now let us assume that the eigenvalues oi, . . . , amn of A{z) are divided into n 
groups of m numbers: 

{ai, . . . , amn] = {a['\ . . . , a«} U • • • U {ai"\ . . . , a^J^)}. 

Lemma 3.3 shows that for a generic A{z) we can construct uniquely defined 
Bi, . . . ,Bn e Mat(m, C) such that for any i = 1, . . . , n, Sp{Bi) = {a^\ . . . , Om } 
and z — Bi is a right divisor of A{z).^ By Lemma 3.1, Bi,...,Bn define A{z) 
uniquely up to a left constant factor, because the eigenvectors of Bi must be eigen- 
vectors of A{z). 

Lemma 3.4. For generic Bi, . . . , Bn G Mat(m, C) with Sp{Bi) = {a^j ^}; there 
exists a unique monic degree n polynomial A{z) = + . . . such that z — Bi 

are its right divisors. The matrix elements of Ai, . . . ,An are rational functions of 
the matrix elements of B\, . . . ,Bn and eigenvalues {af'}- 

Remark 3.5. 1. Later on we will show that, in fact, these rational functions do 

not depend on | aj*' } . 

2. Clearly, the condition of A{z) being monic can be replaced by the condition of 
A{z) having a prescribed nondegenerate highest coefficient Aq. 

Proof. The uniqueness follows from Lemma 3.1. To prove the existence part, we 
use induction on n. For n — \ the claim is obvious. Assume that we have already 
constructed A{z) = z"^^ + Ai . . . such that iJi, . . . , B^-i are its right divisors. 
Let {vi} be the eigenvectors of B„ with eigenvalues {a^^^}. Set Wi = A{a^^^)vi 
and take X S Mat(m,C) such that Xwi = a,'f'\vi for all i = l,...,m. (The 
vectors {wi} are linearly independent generically.) Then A{z) = [z — X)A{z) has 
all needed properties. Indeed, we just need to check that z — Bn is its right divisor 
(the rationality follows from the fact that computing the eigenvectors with known 
eigenvalues is a rational operation). For any i = 1, . . . , m, we have 

[B^ + ArB^-^ + ... + An)vi = ((af V + + ■ ■ ■ + A^) Vi 

= {a'r^-X)A{a^^)vi = {a^^-X)wi = 0. 
Lemma 3.2 concludes the proof. □ 

*It is obvious that the condition on A(z) used in Lemma 3.3, is an open condition. The 
corresponding set is nonempty because it contains diagonal A{z) with {a^*^'} being the roots of 
Akk{z). Similar remarks apply to all appearances of the word "generic" below. 
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Thus, wc have a birational map between matrix polynomials A{z) = Aqz^ + . . . 
with a fixed nondegenerate highest coefficient and fixed mutually distinct eigen- 
values divided into n groups of m numbers each, and sets of right divisors {z — 
Bi, . . . , z — Bn} with Bi having the eigenvalues from the ith group. We will treat 
{Bi} as a different set of coordinates for A{z). 

It turns out that in these coordinates some multipliers R{z) of Theorem 2.1 
take a very simple form. We will redenote by k^'^ the numbers ki, . . . , Kmn used 
in Theorem 2.1 in accordance with our new notation for the eigenvalues of A[z). 
Denote the transformation of Theorem 2.1 with 

n 

K^''=-ki€'L, i = l,...,n, j = l,...,m; Si = ■ ■ ■ = Sm = '^ki 

1=1 

by iS(fci, . . . , fc„). 

(«) 

Proposition 3.6. The multiplier R{z) for 5(0, . . . , 0, 1 , 0, . . . , 0) is equal to the 
right divisor z Bi of A{z) corresponding to the eigenvalues a^*"*, . . . ,am • 

Proof. It is easy to see that if Bi has eigenvalues al ,. .. , am , and z — Bi is a right 
divisor of A{z) then R{z) = z — Bi satisfies all the conditions of Theorem 2.1. 

Conversely, if R{z) is the corresponding multiplier then R{z) is a product of 
n elementary multipliers with one k equal to —1 and one 6 equal to +1. The 
explicit construction of the proof of Theorem 2.1 shows that all these multipliers 
are polynomials, hence, R{z) is a polynomial. The fact that Si = ■ ■ ■ = Sm implies 
that R{z) is a linear polynomial of the form z — B foY some B G Mat(m, C) (to see 
that, it suffices to look at the asymptotics of the canonical solutions). We have 

A{z) = R-^{z + l)A{z)R{z) = {z + I- B)-'^A{z){z - B). 

Comparing the determinants of both side we conclude that Sp{B) = {a^*\ . . . , am}- 
Since no two eigenvalues are different by an integer, B and B — I have no common 
eigenvalues. This implies that (z + 1 — B)~^A{z) must be a polynomial, and hence 

z — _B is a right divisor of A{z). □ 

For any k = (fci, . . . , fc„) G Z" we introduce matrices Bi{k),. . . , Bn{k) such that 
the right divisors of S{ki, . . . , kn)A{z) have the form z — Bi{k) with 

Sp{B,{k)) ^ {a''(^ -h,..., a« - ki\, i^l,...,n. 

They are defined for generic A[z) from the varieties M{---) introduced in the 
previous section. 

Proposition 3.7 (difference Schlesinger equations). The matrices {B,(fc)} 
(whenever exist) satisfy the following equations: 

Bi{...)-Bi{...,kj + l,...) = Bj{...)-Bj{...,ki + l,...), (3.1) 
Bj{. ..,h + l,... )Bi{. ..)= B,{. ..,% + !,... )B,{. . . ), (3.2) 
Bi{ki + l,...,k„ + l) = Ao^Biiki, k„)Ao - I, (3.3) 

where i,j = 1, . . . ,n, and dots in the arguments mean that other ki 's remain un- 
changed. 
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Remark 3.8. The first two equations above are equivalent to 



{z - Bi{. ..,kj + l,...)){z-Bj{...)) = {z-Bj{...,ki + l,...)){z-Bi{...)). 

(3.4) 

Proof of Proposition 3.7. The uniqueness part of Theorem 2.1 imphes that 

{i) 0) 

5(0,...,0, l,0,...,0)o5(0,...,0, l,0,...,0)o5(fci,...,fc„) 

= 5(0, . . . , 0, 1 , 0, . . . , 0) o S(0, . . . , 0, 1 , 0, . . . , 0) o S{ki , . . . , /c„). 

Thus, the corresponding products of the multipliers are equal, which gives (3.4). 
This proves (3.1), (3.2). The relation (3.3) follows from the fact that the multiplier 
for 5(1, . . . , 1) is equal to Aq^A{z), imd A{z) = 5(1, . . . , 1) = A^^A{z + l)Aa. 
This means that the right divisors for A{z) can be obtained from those for A{z) by 
shifting z by 1 and conjugating hy Aq. □ 

Theorem 3.9. Fix ran complex numbers such that no two of them 

are different hy an integer, and an integer M > 0. Then for generic Bi, . . . ,Bn G 
Mat(m, C), Sp{Bi) = {0^*^}™!^, there exists a unique solution 

{Bi{ki, . . . ,kn) : max \ki\ < M} 

i—l,...,n 

of the difference Schlesinger equations (3.1)-(3.3) with 

Sp{Ao) = {pi,..., pn}, Pi/Pj ^ K for ij^j, Pi for i = l,...,n, 
such that 

Sp{Bi{ki, . . . , kn)) = Sp{Bi) — ki and Bi{0, . . . ,0) = Bi for all i = 1, . . . ,n. 

The matrix elements of Bi{k) are rational functions of the matrix elements of the 
initial conditions {5^}"^^. Moreover, these rational functions do not depend on the 

eigenvalues lap-*}. 

Remark 3.10. As we will see later, this theorem also extends to the case of arbi- 
trary invertible ^o- 

Proof. The existence and rationality of the flows have already been proved. Indeed, 
without loss of generality we can assume that is diagonal (the equations (3.1)- 
(3.3) remain intact if we conjugate all Bi{k) and Aq by the same constant matrix). 
By Lemma 3.4 we can construct a (unique) degree n polynomial A{z) with the 
highest coefficient Aq, such that {Bi} is the set of its right divisors. Then, using 
Theorem 2.1, we can define S{k) and hence {Bi{k)}. By Proposition 3.7 they 
will satisfy (3.1)-(3.3). Moreover, all operations involved in this construction are 
rational. 

Thus, it remains to prove uniqueness and the fact that the rational functions 
involved do not depend on the eigenvalues. A simple computation shows that for 
any X, Y,S,T€ Mat(m, C), the relation {z -X){z-Y) = {z- S){z - T) implies 

Y = {X-S)-'^S{X-S), T = {X - Sy^XiX - S), (3.5) 

X = {Y -T)T{Y -T)-\ S = {Y -T)Y{Y -T)-\ (3.6) 
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whenever the corresponding matrices are invertible, cf. [GRW]. Applying this ob- 
servation to (3.4), we see that, generically, {Bi = i?j(0, . . . , 0)} uniquely define 
all 

B,(e«,...,e« ,0,e«„...,e«), f = 0, 1. (3.7) 

Moreover, they are all given by rational expressions involving the initial conditions 
{Bi} only. To move further, we need the following lemma. 

Lemma 3.11. For generic X,Y e Mat(TO, C) with fixed disjoint spectra, there 

exist unique 5, T € Mat(m, C) such that 

{z-X){z-Y) = {z-S){z-T), Sp{S) = Sp{Y), SpiT) = SpiX). 

The m,atrix elements of S and T are rational functions of the matrix elements of 
X and Y which do not depend on the spectra of X and Y . 

Proof 1. Lemma 3.3 proves the uniqueness and shows how to construct T is we know 
the eigenvalues x\,. . . ,Xm oi X and vectors Vi such that (xj — X){xi — Y)vi = 0. If 
wc normalize ViS in the same way, for example, by requiring the first coordinate to 
be equal to 1 (this can be done generically), then using the construction of Lemma 
3.3 we obtain the matrix elements of T as rational functions in the matrix elements 
of X, Y and .xi, . . . , a;„. However, it is easy to sec that these rational functions are 
symmetric with respect to the permutations of xi, . . . , j:„, which means that they 
depend only on the elementary symmetric functions X]ii<-<ifc ^ii"' ^ik of 
But these are the coefficients of the characteristic polynomial of X, and hence they 
are expressible as polynomials in the matrix elements of X. □ 

Proof 2, see [O]. The uniqueness follows from Lemma 3.3. To prove the existence, 
denote by A the solution of the equation FA — AX = I. Generically, it exists, it is 
unique and invertible. Set 

S = X + A-\ T = Y-A-\ 

Then it is easy to see that {z — X){z — Y) = {z — S){z — T). Furthermore, if 

Y and T have a common eigenvalue then they must have a common eigenvector, 
which contradicts the invertibility of F — T = A^^. Hence, Sp(T) = Sp{X) and 
Sp{S) = Sp{Y). □ 

Remark 3.12. In the case of 2 by 2 matrices, it is not hard to produce an explicit 
formula for S and T in terms of X and Y: 

S ={X + Y-Tr:Y)Y{X + Y-TiY)-'^, T = {X + Y -Ti X)-'^X{X + Y -Tr X). 

(3.8) 

Now let us return to the proof of Theorem 3.9. Recall that we already proved 
that the initial conditions define (3.7) uniquely. Now let us use (3.4) with 

0) 

(fci, . . . , A;„) = (1, . . . , 1, , 1, . . . , 1), j^i. 

By (3.3), we know what Sj(l, . . . , 1) is. Thus, we know both matrices on the left- 
hand side of (3.4), and hence, by Lemma 3.11, we can compute both matrices on 

U) 

the right-hand side of (3.4), in particular, Bj(l, . . . , 1, , 1, . . . , 1). 
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Now take (3.4) with 

(i) (0 
(fci,...,fc„) = (l,...,l, 0,1,...,1, 0, !,...,!), 

where i,j,l are pairwise distinct. Applying Lemma 3.11 again, we find all 

0) (0 
i3,:(l,...,l, 0,L...,1, 0,1,...,1). 

Continuing the computations in this fashion (changing one more 1 to in (A:i, . . . , fc„) 
on each step), we obtain all 

iiHei ,ei_i, 1 ••• -0,1. 

Together with (3.3) (and (3.7)) this computes all Bi{k) with max \ki\ < 1. Iterating 
this procedure, we complete the proof. □ 

4. An alternative description of the Schlesinger flows. 

The goal of this section is to provide yet another set of coordinates for the 
polynomials A{z), in which the flows described in the previous section can be easily 
defined. In particular, this will lead to a different proof of Theorem 3.9, which will 

be valid for an arbitrary invcrtible Aq. 

Proposition 4.1. In the assumptions of Theorem 3.9, the monic degree n polyno- 
mial 

[z - Bi(0, 1, . . . , l))(z - i?2(0, 0, 1, . . . , 1)) • • • (^ - B„(0, . . . , 0)) 

has z — Bi, i = 1, . . . ,n, as its right divisors. 

This statement and Theorem 3.9 provide a proof for Remark 3.5(1). 

Proof. Using (3.4) we obtain (j > i) 

U) \ ( (i+i) 

^-Bi(0,...,0, l,...,l)j f^-B,(0,...,0, 1 ,...,1) 

(i) + 1) \ / + 1) \ 

z - 5,(0,..., 0, 1,0,..., 0, 1 ,...,l)j (z-B,(0,...,0, 1 ,...,l)j. 

Using this commutation relation, we can move the factor (z—i?j (•••)) in the product 
above, to the right most position, where it will turn into [z — 5^(0, . . . , 0)) = 
{z-Bi). □ 

Let us introduce the notation (Zi , . . . , Z„ e Z) 

Ci(Zi, ...,/„):= Bi(Zi, . . . + 1, . . . , Z„ + 1), i = 1, . . . ,n. 

If we denote by A{z) the polynomial of degree n with highest cocfficic^nt Aq such 
that {z — Bi} are its left divisors, then the definition of Bi(k) and Proposition 4.1 
imply that for any / = (li, . . . ,ln) G Z", 

S{h, . . . , ln)A{z) =A^{z- Ci{l)) ■■■{z- Cn{l)). (4.1) 

(To apply Proposition 4.1, we also used an easy fact that for any solution {Bi{k)} 
of (3.1)-(3.3) and any Zi, . . . , Z„ e Z, 

B[{ki, . . . ,kn) := Bi{ki + Zi, . . . , fc„ + Z„), i = l,...,n, 

also form a solution of (3.1)-(3.3).) 
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Lemma 4.2. The map {Bi} i— > {Cj} is birational. 

Proof. The rationality of the forward map follows from Theorem 3.9. The ratio- 
nality of the inverse map follows from Lemma 3.3 (indeed, we just need to find the 
right divisors of the known matrix S{li, . . . , ln)A{z)). Even though it looks like to 
construct Bi wc need to know the eigenvalues of Ci , it is clear that by normalizing 
the eigenvectors of A{z) corresponding to these eigenvalues, in the same way, we will 
obtain a formula for Bi which will be symmetric with respect to the permutations 
of these eigenvalues, and thus we can rewrite it through the matrix elements of Cj's 
only (this argument was already used in the first proof of Lemma 3.11 above). □ 

Our goal is to describe the transformations S{k) in terms of {Ci}. We need a 
preliminary lemma which generalizes Lemma 3.11. 

Lemma 4.3. For generic Xi,. . . , Xn e Mat(m, C) with fixed disjoint spectra and 
any permutation a € Sn, there exist unique li,...,YAr € Mat(m, C) such that 
Sp{Yi) = Sp{Xi) for alli^l,...,N, and 

_ Xi) • • • (^ - Xjv) = {z- ■■■{z- y<,(jv)). 

The matrix elements of {Yi} are rational functions of the matrix elements of {Xi} 

which do not depend on the spectra of {Xi}. 

Proof. The existence and rationality claims follow from Lemma 3.11, because ele- 
mentary transpositions {i,i+l) generate the symmetric group ^jv- To show unique- 
ness, we rewrite the equality 

_ YD ...[z-Yl,) = iz- F/O ■■■{z- Yi^), Sp{Yl) = Sp{Yl'), 

in the form 

{z - Yir\z - Yi) = [{z - y,") ---{z- F^)) ((z - Yi) . . . (z - r;,))"'. 

If the spectrum of F/' is disjoint with the spectra of ^2 7 • ■ • j ^rni then both sides of 
the last equality are entire because they cannot possibly have common poles. Since 
both sides tend to / as 2; — > 00, by Liouville's theorem we conclude that both sides 
are identically equal to /, and F/ = F/'. Induction on N concludes the proof. □ 

Proposition 4.4. In the assumptions of Theorem 3.9, {Ci{l)} satisfy the equations 

{z + l-Ci)---{z+l-Cn)Ao{z-C{)---{z- 
= {z + l-Ci+i)---{z + l-Cn)Ao{z-Ci)---{z-Ci), (4.2) 
Q = • • • I ^n), Cj = Cj{li, . . . , li-i,li + 1, . . . ,ln) for all j, 

and 

Ci{h + l,...,ln + l) = Ao^Ciih, ln)Ao - I. (4.3) 
In both equations i= 1, . . . , n is arbitrary. 

Proof. The relation (4.3) is a direct corollary of (3.3). Let us prove (4.2). Proposi- 
tion 3.6 implies that the multiplier for the shift 



{h, ■ ■ ■ ,ln) '-^ {h, ■ ■ ■ , h-lji + 1, h+l, ■ ■ ■ ,ln) 
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has the form R{z) = z — Bi{l). Thus, (4.1) gives 



Aa{z - Ci) • • • - C^){z - BS)) = (z + 1 - BS))Ao{z -C^)■■■{z- C„). 
Comparing the spectra of factors on both sides and applying Lemma 4.3, we get 
Aq{z -C^)---{z-Ci) = {z + l- Bi{l))Ao{z -Ci)---{z- Ci_i), 

{Z - Ci+i) ■■■{Z- Cn){z - Bi{l)) = {Z-Ci)---{Z- Cn). 

Combining these two relations and shifting z i-^- z + 1 in the second one, we arrive 
at (4.2). □ 

Theorem 4.5. Fixmn complex numbers {aj*''}"l7j=i '^^^'^ ^^^^ '^^ of them are 

different by an integer, an integer M > 0, and any nondegenerate Aq € Mat(m, C). 
Then for generic Ci,...,C„ € Mat(m, C), Sp{Ci) = {aj'^}^_^, there exists a 
unique solution 

{Ci{h,...,ln) ■■ max < M} 

i— l,...,n 

of the equations (4.2) and, consequently, (4.3), such that 

Sp{Ci{h, . . . ,ln)) = Sp{Ci) - k and Ci(0, . . . ,0) = Cj for all i = 1, . . . ,n. 

The matrix elements of Ci{l) are rational functions of the matrix elements of the 
initial conditions {Cij^^, and these rational functions do not depend on the eigen- 
values {ap^}- Moreover, 

Bi{ki, . . . , kn) ~ Ci{ki, . • . , ki, ki-\-i 1, . . . , kn 1), i = 1, . . . , n, 

solve the difference Schlesinger equations (3.1)-(3.3). 

Remark 4.6. If the ratios of eigenvalues of Aq are not real then Theorem 3.9, 
Lemma 4.2, and Proposition 4.4 provide a proof of Theorem 4.5. However, our goal 
is to provide an independent proof of this theorem, thus giving a different proof of 
Theorem 3.9 with arbitrary invcrtiblc ^o, cf. Remark 3.10. 

To prove Theorem 4.5 we will develop a rather general formalism. 

(a) Semigroup. Let P be a semigroup and Po be its subset. We assume that 
every element of Pq has a type. The types of two different elements pi,P2 G Po may 
be the same, which will denoted by t(pi) — t{p2), and may be disjoint, which will 
be denoted by -L t{p2). The types may also be neither equal nor disjoint. 

Assumption 4.7. For any elements pi,...,pN G Pq such that their types are 
pairwise disjoint: 

t{pi)±t{pj), ij^j, l<i,j<N, 

and for any permutation a E Sn there exist unique elements pi, ■ ■ - Pn € -Po such 
that t{pi) = t{pi), i = 1, . . . ,N, and 

P1---PN =Pail)---PaiN)- 
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We will be interested in the situation when 

p ^ pMat(m,C) ^ l^fe ^ Q^^k-l + . . . + I g. g Mat(m, C), 1 < i < fc} , 
Po = pMat(m,C) ^ |^ _ g | g g Mat(m, C)} , 

t{z-Q) = {z&C\ Aet{z - Q) = 0} = ^^(Q). 

The notions of equality and disjointness for types are the natural ones for the 
m-point subsets of C. Lemma 4.3 shows that pMat(m,c) pMat(m,c) ^^^j^j^ 
Assumption 4.7 generically. 

Proposition 4.8. Let P he a semigroup satisfying the Assumption J^.l. Assume 
that we have an equality in P of the form 

where all p^^ \qf'^ are from Pq, the types of all elements on the left hand side 
are pairwise disjoint, the types of all elements on the right-hand side are pairwise 
disjoint, and 

{i(p(^)),...,i(p«)} = {i(«?'),---,<(aH])} 

for all j = 1, . . . , k. Then 

Proof. By Assumption 4.7, for any j = \, . . . ,k, there exist . . . , qm] such that 

and t{pf^) = t{S[f '^). Then by the uniqueness part of Assumption 4.7 applied to 
(4.4) we obtain p^p — qf ^ for all i,j. This immediately implies the claim. □ 

We, essentially, used Proposition 4.8 in the proof of Proposition 4.4 above. 

(b) Commuting flows on sequences. Denote by V the set of all sequences 
{pk}kez C Po such that the types of all elements of a sequence are pairwise disjoint. 

Fix an integer integer n > 0. For any I G Z we define a map Fi : P ^ V as 
follows: 

Fi : {pk}kez {Qk}kez, 

Pl+fmPl+fin+l ■ ■ ■ Pl+(fi+l)n-l = Ql+iJ.n+l(}l+fm+2 ' ' ' Ql+{fj.+l)nj M S Z, 

t{pj-n), if (j - divides n, 
t{pj), otherwise. 



In this definition /z ranges over all integers, and for each fj, we use Assumption 4.7 
for (J = {12 ■ ■ ■ n) £ Sn- Clearly, Fi is invertible. 
It is convenient to denote 

Fh{Fi,{---Fi^ {{Pk}))) = {p'^'-''-}. 
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Then the second line in the definition above takes the form 

Pl+fj.nPl+p,n+l ■ • • Pl+{fj,+ l)n-l = p'l+fj.n+lP'l+tJ.n+2 ' ' ' p\+{^+l)n- 

For example, for n = 2 we have 

P2s-lP2s = pIsPL+I^ t{P2s) = t{P2s), t{pls+l) = t{p2s-l), S £ Z, 

P2sP2s+l = pIs+iP2s+2, i(pL+l) = *(P2s+l), i(pL+2) = *(P2s), S € Z. 

It is immediately seen from the definition that Fi+^n = Fi and 

1,1 + 1,. ..,l+n-l (. r\ 

Pk = Pk-n (4.5) 

for any fc, Z G Z. 

Theorem 4.9. (i) For any i,j £ Z, F^ and Fj commute. That is 

Pk^ = Pk' f'^^ ""2/ {Pk} e v. 

(ii) For any {g^} G V and any i,j £ Z such that < j — i < n, set pk = q]^^''"''' . 
Then 

(lipj=pUi- (4-6) 

Remark 4.10. Part (i) of this theorem means that we have defined an action of 
Z" on v. There is a much larger group which acts on P. Let tt : Z ^ Z be a 
bijection such that for any fc G Z the sets 

Ik = {i e Z : i < k, 7r(i) > 7r{k)}, Jfc = { j e Z : j > k, 7r(j) < 7r(A;)} 

are finite: / = . . . ,is}, J = {ji, . ■ ■ ,jt}- Then, given a sequence {pk} S V, we 
define {p^} G P by 

Ph ■■■Pis Pk Ph ■ ■■Pit = P'h ■ --Ph Plik) P'h ■■■Pis 

where t{p'i) = t{pi) and = t{pk). One can show that this defines an action 

of the group of all n satisfying the condition above on the space V. The maps {Fi} 
correspond to shifts by n along n nonintersecting arithmetic progressions {I + fin : 

/i G Z}, hence they must commute. 

Proof of Theorem 4-9. Since Fi = -F(+„, it suffices to assume that < j — i < n. 
Consider the product 

^=PiPi+l---Pj+2n-l- 

On one hand, we have 

11= {pi- ■■Pi+n-l) {pi+n ■ ■ ■Pi+2n-l)Pi+2n ■ ■ ■ Pj+2n-l 
= {pI+i ■ ■ -pi+n) {pl+n+l ■ ■■pi+2n)Pi+2n ■ ■ ■ Pj+2n-l 
= Pi+i ■ ■ ■ (p* • • • Pj+n-l)P)+n ■ ■ ■ P\+2n Pi+2n ■ ■ ■ Pj+2n-l 

= Pi+i ■ ■ ■ pli'li ■ ■ ■ P^'ln pi+n ■ ■ ■ Pl+2n Pi+2n ■ ■ ■ Pj+2n-l- 
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On the other hand, we have 

= Pi--- Pj-l • • • {plj+n+l ■ ■ ■ P'j+2n) 

= Pi - ■ -Pj-l p'j+l ■ ■ -Pi+n-l {Pt+n ' ' 'Pi+2n-l)pi+2n ' ' ' Pj+2n 

j j i^j i^j j i 

— Pi-- -pj-\Pj_^i - - - P^j^n-lPi+n+l ' ' ' Pi+2nPi+2n ' ' ' Pj+2n- 

Thus, we obtain 

{pI+1 ■ ■ -P)-! • • -P^iln) • • -P^jin) {p)+n ' ' ' P\+2nPi+2n ' ■ ■ Pj+2n-l) 

~ \Pi' ' 'Pj-lPj^l ' ' ' Pi^n-l) [Pi^n-^l ' ' ' Pj+n ) \Pj-\-7i-\-l ' ' ' Pi-\-2n Pi-\-2n ' ' 'Pj-\-2n)- 

(4.7) 

Comparing the types in the three factors on the left and on the right, we see that 
we are in a position to apply Proposition 4.8. It implies, in particular, that the 
middle factors arc equal. Since the order of the types in the middle factors is the 
same, these middle factors must be equal termwise: 

p'k'=t^k^ i + n+l<k<j + n. 

Since Fi = Fi+n for all n, and i and j are arbitrary, we see that p^.'* = p]l-' , i + 1 < 
k < j. Switching from (i, j) to {j, i + n), we get p^j-' = p^'* for j + l<k<i + n. 
Thus, the commutativity relation is proved for i + l<k<i + n, and thus for all 
k G Z. The proof of the first part of Theorem 4.9 is complete. 

In order to prove Theorem 4.9 (ii), we need to compare the first and the third 
factors of the two sides of (4.7). The first factors give 

i i j^i j.i j j 

Pi+1 ■ ' 'Pj-lPj + 1 ■ ■ 'Pi+n — Pi' ' 'Pj-lPj+l • • -Pi+n-l- 

Commuting pi in the right-hand side to the right and using Proposition 4.8, we see 
that 

Pi---Ps=P\+i---p\Xs,i (4.8) 

where i + 1 < s < j — 1, Xs,i G Pq, and t{xs^i) ~ t{pi). Since j is arbitrary (but 
j — i < n), we can assume that i + 1 < s < i + n—2. Note also that (4.8) also holds 
for s = i + n — 1 with Xi+n-i,i = P\+n^ as follows from the definition of Fj. 

Similarly, looking at the third factors and substituting {pk} by F^^{pk}, we get 

Pt ■ --Pj = VL, {Fr^iPk}), ■ ■ ■ {Fi'{Pk})^_, (4.9) 

where j — n + 2<t<j — 1, ytj G Pq, and t{ytj) = t{pj). Again, this also holds 
for t = j - n + 1 with yj+„_ij = {F~^ {pk}) 

Lemma 4.11. For i + l<s<i + n— 1, we have 

Xs,i = {F-^\o---oF-'{pk}).. 

Proof. Induction on s — i. To prove both the base s = i + 1 of the induction and 
the induction step we first use (4.9) to write 

Pi---Ps = yi,s {F-HPk}),- ■ ■ {Fs~HPk}),_, 
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and now use the induction hypothesis on the factors which go after ?/j_s to obtain 

P^ ■ ■ -Ps = yM {F^F-\pk})^^, ■ ■ ■ {F^F-'{pk})^_, {Fr^\ o • • • o F-\pk}). . 

(If s = i + 1 then the second step is empty.) Since 

t {{Fr^\ o • • • o F-^{pk}).) = t{xs,i) = t{pi), 

comparing with (4.8) we conclude that Xg^i = {F~^-^ o • • • o F~^{pk}).. 

This argument works for s < i + n — 2. For s = i — n + 1 the lemma follows from 
(4.5). □ 

Now we return to the second part of Theorem 4.9. Applying Lemma 4.11 to all 
but one factors in pi ■■ ■ pj , and then to all factors in • • • pj , we obtain 

Pi---Pj= pi+i ■ • -p}.! {Fri\ o • • • o Fr_\{pk}).Pj 
= pi+,...p^_,p) {F-+\o...oFr^{pk}).. 

In the last two products all but last two factors coincide. By Proposition 4.8, this 
means that the products of last two ones also coincide: 

{Fr^\ o • • • o Fr_\{p^})^p^ = p\ {Fr^\ o • • • o Fr^{pk}). . 

Renaming F^-^ o ■ ■■ o F~^{pk} by {qk} we arrive at (4.6). The proof of Theorem 
4.9 is complete. □ 

(c) Proof of Theorem 4.5. Let us concentrate on the case P = p^iat(™.c)^ 
(a) above. Since Assumption 4.7 generically holds in this case 
(see Lemma 4.3), we will be acting like it always holds, keeping in mind that all 
the claims we prove hold only generically. Set Pi = z — Ci, i = l,...,n, where 
Ci = C'i(0, . . . , 0) are as in Theorem 4.5. More generally, define 

p^+^,n^z~fi-A^aiO,...,0)AQ^, i^l,...,n, ^ e Z, (4.10) 

where Aq is an arbitrary invertiblc clement of Mat(m, C). The assumption that no 
two numbers of the set {aj'''} are different by an integer, guarantees that {pk} € 'P- 
Now define {Ci{h, . . . , In)} by 

(Fl'---Fi^"{pk}) =z-n-A^Ci{h,...,ln)Ao'', ^=l,...,n, /x e Z. 

(It is immediately seen that the subset of V consisting of sequences {pk = z — Qk} 
such that Qk+n = I + AoQkAg^ is stable under the flows Pi, ... , Fn, which shows 
that Ci{l) are well-defined.) The very definition of Fj implies (4.2). Furthermore, 
(4.3) is a direct corollary of (4.5). It is easy that Sp{Ci{l)) = Sp{Ci) — k, and 
this gives the existence part of Theorem 4.5. The uniqueness and rationality claims 
follow from Lemma 4.3. Finally, let us show that 

Bi{h\, . . . , hn) = Ciikli . . • , ki, ki^\ 1, . . . , kn 1) 
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solve (3.1)-(3.3). The relation (3.3) is equivalent to (4.3). We will derive (3.4) (and 
hence (3.1), (3.2)) from Theorem 4.9(ii). 
Fix 1 < i < j < n and define 



Then 





{Pk}=F^^- 


■■Fn-iPk}, 


{»} = 


F-^\o...oF-'{pk}. 


Pj 


= z- Cj{ki, . . 


■ 7 kn) = Z 


Bjiki,.. 


■ ,kj, kj-^-i + 1,...,A;„ + 1), 




= z- Bi{ki, . . 


■ 1 kj, + 


1 k 


+ 1)7 




= z - Bj{ki, . . 


■ 7 ki—ij ki + 


1, ki+i, . . 


■ , kj , kj^i + 1, . . . , fc„ + 1). 




= z- Bi{ki, . . 


• ' ^j—i ' ~^ 


1, . . . , kfi 


+ !)• 



If we apply the shift 

kj-\-l I ^ ^i+l I7 • * • 7 ' ^ kft 1, 

then the equality g^p^ = p*gj turns into (3.4). This completes the proof of Theorem 
4.5. □ 

Remark 4.12. The set of sequences {pk — z — Qk} with Qk+n — ^ + -^oQk-^o 
is also stable under the action of permutations it I'E^ Z (see Remark 4.10) of the 
form 

n{i + fin) = a{i) + fin, a G Sn, i = 1, . . . ,n, fj, gZ. 

Define Ci{l) = Ci{l) + lil. Then we obtain a birational action of the semidirect 
product Z" K Sn on {C\, . . . S (Mat(m, C))" which preserves the spectra of 
C'iS. 

Remark 4.13. If instead of (4.10) we use periodic initial conditions 
Pi+nn = z-Ci, i = l,...,n, fiGZ, 

which corresponds to the autonomous limit of the difference Schlesinger mentioned 

in the Introduction, then the maps Fi, . . . , Fn are exactly the monodromy maps 
constructed by Veselov [V] in the framework of set-theoretical solutions of the 
quantum Yang-Baxter equation. We refer to [V] for details and further references 
on the subject. 

Remark 4.14. Solutions of the g-difference Schlesinger equations mentioned in the 
Introduction are obtained from considering {pk = z — Qk} with Qk+n = Q^oQk^o^ ■ 

5. Continuous limit 

We start with a brief survey of the classical deformation theory for linear matrix 
differential equations, which is due to Riemann, Schlesinger, Fuchs, and Garnier, 
see [JMU], [JM] for details. 

Consider a first order matrix system of ordinary linear differential equations 

^=B(C)y(C), BiO^B^ + j^TZ^- (5.1) 
at, c, Xk 
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Here all matrices are in Mat(m, C). We will assume that all B^s can be diago- 
nalized: 



Bk = GkTkGl\ Tk = diag(4^\ . . . , i^), k = l,...,n, 
Boo = Goo diag(si, . . . , Sn)G^ 

where t'^''^ — ^"L, i ^ j, for all k ^ oo, and Sj ^ s^-, i ^ j. 
Alternatively, we may also consider 

in which case we assume (in addition to the above assumption on S^, A; = 1, . . . , m) 
that 

n 

- = GooTooG-\ Too = diag(4°°\ . . . 

fe=i 

with t^°°^ - 4°°' ^ Z for i ^ j. 

Since we can conjugate y and {Sfc} by Gqo, we may set Goo = i without loss of 
generality. 

One can show, see e.g. [JMU, Proposition 2.1], that there exists a unique formal 
solution y{Q of (5.1) or (5.2) of the form 



3^(C) = >'(c)cxp(r(c)), >'(c) = /+i>ir'+i>2r' + ---, (5.3) 



where 



r(C) = 



diag(si, . . . ,s„)^ + Tex, ln(^; ^) for (5.1), 
Tooln(^-i) for (5.2), 



with Too = diag(ti°°\ . . . ,t^'). (This formula is also the definition of Too for 
(5.1)). This is the analog of Propositions 1.1 and 1.2. 

It turns out that for (5.2) the series in (5.3) is convergent, and after multiplication 
by exp(T(C)) it defines a holomorphic near ^ = oo multi-valued function y°°{Q. 
However, in the case of (5.1) this series is, generally speaking, divergent. Then the 
analog of Theorem 1.3 holds. Namely, there exist unique holomorphic solutions 
y^'"^ of (5.1), defined for 3?^ <C and ^ 0, respectively, such that they have 
the asymptotic expansion (5.3) as C ^ co-^ Since both these functions solve the 
same differential equation, there exist constant matrices such that the analytic 
continuations of J^''*" in 3C » (-C 0) are related by y^ = y^S"^. The matrices 
are called the Stokes multipliers. 

It is also possible to determine the nature of solutions of (5.1), (5.2) near the poles 
Xi,. . . ,Xk- Namely, one can show that there exist locally holomorphic functions 

y^'^ {0 = 1+ yi'^ {z-x,)+ yi'^ {z-xkf + ... 



^As in the case of difference equations, one has to be careful in choosing the sector where 
( may tend to oo. One may always taJje argC^ € {it/2 + e,37r/2 — e) for SRC <S and argf G 
(— 7r/2 + e, 7r/2 — e) for 3?^ 2> 0. If Si — sj ^ M, these sectors may be extended. 
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such that for any solution y{C), there exist constant matrices Ck such that locaUy 
near = Xk we have 

y{0 = Gfej)W(C)exp(Tfcln(^ - Xk))Ck. 

(Recall that {Gk} were defined above by Bk = GkTkG^^.) In particular, if we fix 
paths from ^ = oo (or ±00 for (5.1)), then we can define {Cfe} for the (analytic 
continuations of the) canonical solutions 3^°° or J^''*". 

Thus, to any equation of the form (5.1) or (5.2), we associate the following 
monodromy data: {Tk}^^i and Too, {Ck}k=i computed for the canonical solution 
y°° or y^'^ , and in the case of (5.1) we also add the Stokes multipliers 5^ and the 
exponents si, . . . , Sm. 

If we analytically continue any solution y{() of (5.1) or (5.2) along a closed path 
7 in C avoiding the singular points {xk} then the columns of y will change into 
their linear combinations: y 1-^ yMj. Here is a constant invertible matrix 
which depends only on the homotopy class of 7. It is called the monodromy matrix 
corresponding to 7. If 7 is a positive loop around Xk then the corresponding mon- 
odromy matrix Mk for the canonical solution 3^°° or 3^'''" can be computed using 
the monodromy data introduced above: 

Mk = G^^ cxp(2mTk)Ck. 

The basic problem of the isomonodromy deformation of the linear system (5.1) or 
(5.2) is to change B{() in such a way that the monodromy data, or, more generally, 
the monodromy matrices {Mk} remain invariant. 

There are two types of isomonodromy deformations, both were discovered by 
Schlesinger [Sch] and later generalized to the case of singularities of higher order in 
[JMU], [JM]. 

The first type is a continuous deformation which allows the singularities xi, . . . ,x„ 
to move and describes {Bk} as functions of xj's. This deformation leaves the whole 
set of monodromy data intact. The evolution of {Bk} is described by a system of 
partial differential equations called the Schlesinger equations: 

dx-^x -xj' db- ~ ^ xj-x- ^"-""^J' J>« - ■^'•••''^> 10-4; 

l<l<n 

where for (5.2) the term with B^ is dropped. 

It is not hard to show that this system has a local solution for arbitrary initial 
conditions {Bk{x\, . . . ,x°^)}. It is a much deeper fact (proved independently in 
[Mai], [Miw]) that the Schlesinger equations with arbitrary initial conditions have 
a global meromorphic solution on the universal covering space of 



{(a;i, . . . ,Xn) & : Xi ^ xj for i ^ j}. 



To describe this fact, one often says that the system of Schlesinger equations enjoys 
the Painleve property. 

The second deformation (or, better to say, transformation) is an action of 
on the space of B{z), which consists of multiplying y{z) by an appropriate rational 
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function on the left: y{z) ^ Tl{z)y{z). Such a transformation (called Schlesinger 
transformation) is uniquely determined by the shifts 

tf^tf^+xf\ k = l,...,n,oo, 

of the eigenvalues of Bk- Here all A^*^^ are integers, and their total sum is equal to 
zero. Schlesinger transformations exist for generic {Bk}, see [JM]. Clearly, these 
transformations change the monodromy data, but they do not change the mon- 
odromy matrices {M^} and the Stokes multipliers S^. 

Now let us take a difference equation of the type considered earlier: 

Y{z + l) = A{z)Y{z), A{z) = Aoz"" + Aiz"-' ■ ■ ■ + An. (5.5) 

We distinguish two cases (cf. Propositions 1.1 and 1.2): 

• Aq is diagonal and has pairwise distinct nonzero eigenvalues; 

• Aq ~ I, Ai is diagonal and no two eigenvalues of Ai are different by an integer. 

As was explained in §4, see Proposition 4.1, we can generically represent A{z) 
in the form 

A{z) = Ao{z-C\)---(z-Cn), 

where the eigenvalues {a^*'} of {Cj} are zeros of det A{z) divided into n groups 
of m numbers. We assume, as usual, that no two eigenvalues are different by an 
integer. 

Suppose that A{z) depends on a small parameter e, and as e ^ 0, we have 
Ci - Vie-'' + Bi ^ 0, i = l,...,n, {Aq - I)e-^ - B^^ ^ 0, e ^ 0, (5.6) 

for some pairwise distinct complex numbers y\,...,yn and some Bi, . . . , B„, Boo € 

Mat(m,C). (The limit relation for is omitted in the case = I ■) 

Note that if we multiply the unknown function Y{z) in (5.5) by Hi r(z — j/iC"^), 
then (5.5) takes the form 

r(^ + 1 - yie-^)Y{z + 1) j = (7 + B^e + o(e)) 

K-f^)--(-f^)(n^<---H- 

If we now assume that ]"[,• r((C + yi)(r^)Y{C,(r^) tends to a holomorphic hmction 
yiC,), then the difference equation above in the limit e ^ turns into the differential 
equation (5.1) (or (5.2)) with Xi = y,. 

Substituting the asymptotic relations (5.6) into (4.2), we see that for any fixed 
Zi , . . . , ^„ S Z, we have 

Cj(Zi, ...,/„) + - t/je~^ + ^ 0, i = l,...,n, e ^ 0. 

(This conclusion is based on the fact that if X = xe~^ + Xq + o(l), Y = ye~^ + 
Yo + o(l), where x,y e C, x ^ y, and {z - X){z - Y) = {z - S){z - T) with 
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Sp{S) = Sp{Y), SpiT) = Sp{X), then S = Y + o{l), T = X + o{l), see the exphcit 
construction of T in Lemma 3.3.) 
In particular, for any fci, . . . , A;„ € Z 

Bi{ki,...,kn) + ki-yie~''^ +Bi^O, i=l,...,n, e ^ 0. 

(See §4 for the relation of {Bi} and {Cj}.) Thus, on finite intervals Bi{k)+ki—yie~^ 
is approximately constant. However, if wc assume that + fc, — yiC~^} for k 

of size e^^ approach some smooth functions of ekj: 

Bi{[xie~'^], [x„e"^]) + (x, - yi)er^ + Bi{xi, . . . ,a;„) ^ 0, e ^ 0, 
B^{0,...,0)=B^, i = l,...,n, 

then the corresponding equation (5.5) converges to (5.1) with {Bi — Bi{x)} and Xi 
replaced by y, — Xi. Furthermore, the difference Schlesinger equations (3.1)-(3.3) 
tend to 

dB, [Bi,B,] -dBi 

Comparing these equations to (5.4), wc arc led to the following 

Conjecture 5.1. For generic Bi, . . . , BmB^o and pairwise distinct yi, . . . , y„ G C, 
take Bi = Bj(e) e Mat(m, C), i = 1, . . . ,n such that 

Bi{e) - yic-^ +Bi^O, e ^ 0. 

Let Bi(ki, . . . , fc„) be the solution of the difference Schlesinger equations (3.1)-(3.3) 
with the initial conditions {Bi(0) = B^}, and let Bi{xi, . . . ,Xn) be the solution of 
the classical Schlesinger equations (5.4) with the initial conditions 

Bi{yi,...,y„) = Bi, i = l,...,n. 

Then for any xi,. . . ,Xn e M and i = 1, . . . ,n, we have 

Bi([xie"^],.. ., [a;„e"^]) + {x, - ?A)e"^ + B,{yi - xi, . . . , ?y„ - .t„) ^0, e ^ 0. 

As for isomonodromy deformations of the second kind (Schlesinger transforma- 
tions), we are able to prove an asymptotic result rigorously. We will consider the 
case of equation (5.1), for (5.2) the situation is similar. 

Fix 

}i<fe<„,i<,<™CC, t^'^-t^^^Z unless n = j„ k, = k,. 

For any 

BfeCMat(m,C), .5p(^fe) = {tf }™ i, fc = l,...,n, 
and pairwise distinct xi,. . . ,Xn € C, we define 

Bk{e) = Xke~'^ - Bk, k = l,...,n, e ^ 0. 

Then we have 



5p(Bfe(e)) = {af af^=Xke-'-tf\ j = l,. 

We also fix Boo = d.iag(si, . . . , Sm), Si ^ sj for i ^ j. Set 

Ao{e) = I + eBo 
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Lemma 5.2. For generic Bi, . . . , Bn and |e| small enough, there exists a unique 
degree n polynomial A{z, e) = Ao{e)z"- + Ai{e)z"-~^ + ■ ■ ■ having {z — -Bfe(e)} as its 
right divisors. 

Proof. According to Lemma 3.4, the statement is true for large |e|. On the other 

hand, for fixed {Bk} the existence of A{z) is an open condition on e, and if it holds 
for large |e|, it also holds for |e| small enough. □ 

Theorem 5.3. Fix any integers {A^*^ }^i' '^ = ^,■■■■,'^,00, of total sum 0: 

m / n \ 

j=i \i=i / 

Then for generic Bi, . . .Bn and small enough \e\, there exists the transformation of 
Theorem 2.1 for the equation Y{z + 1) = A{z, e)Y{z) with 

K« = -Af, l<k<n, Sj=-X'j°°\ l<j<m. 

Furthermore, if we denote by {Bk} the coefficients of {5.1) after the corresponding 
Schlesinger transformation, and by {_Bfe(e)} the matrices such that {z — Bk(e)\ are 
the right divisors of the transformed A{z,t), then 

ABi(e)A-i -a;ie-i +Bi ^ 0, e ^ 0, z=l,...,n, 

where A = diag(e''^ , . . . , e*" ) . 

Remark 5.4. It is easy to verify the statement of Theorem 5.3 if A^*^ = — Aj°°^ = 
±1 for some fixed i and all j = 1, . . . , m, with all other A's being zero. Then -Bj(e) = 

Bi{e; 0, . . . , 0, ± 1 , 0, . . . , 0). As was mentioned above, for any fixed ki,. . . ,kn S Z, 
we have the asymptotics 

Bi{e;ki, . . . ,kn) + h - Xie~^ + Bi ^ 0, e ^ 0. 

Hence, Theorem 5.3 implies that 

~ _ r Bfc, k^i, 
''~\ Bi±I, k = i. 

This is immediately verified using the fact that the multiplier TZ{C,) for the (contin- 
uous) Schlesinger transformation in this case equals TZ{Q = (C — Xi)^^. 

Proof of Theorem 5.3. Arguing as in the proof of Lemma 5.2, we can show that all 
the statements used in this proof which hold gcnerically (like the existence of the 
polynomial with given right divisors) , also hold for generic Si , .... 6^ and small 
enough e. Thus, we will ignore the questions of genericity from now on. 

Note that we can decompose the transformations of Theorem 5.3 in both discrete 
and continuous cases into compositions of elementary ones of the same type (those, 
for which exactly one of {A^*^} is equal to ±1, and exactly one of {A^^^} is equal 
to =f1) with all others being zero). It is clear that the claim of the theorem follows 
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from a slightly more general claim for the elementary transformations: we assume 
that 

AoBi{e)Ao^ -xie-^ +Bi^O, e ^ 0, l = l,...,n, (5.7) 

with some diagonal Aq containing integral powers of e on the diagonal, and we need 
to conclude that 

AAoBi{e)AQ^A-^ -xie-^ +Bi^O, e^O, l = l,...,n. (5.8) 

Let us consider the elementary transformations with A^-^'' = 1, A^""'' = —1. 
Denote by R{z, e) and 7^(C) the corresponding multipliers for the discrete and con- 
tinuous equations. According to the proof of Theorem 2.1, 

R{z, e) = {z- xie-i + tf^)Ei + Ro{e), R-\z, e) = I - Ei + ^^'^'^ 

z - xie-'^ + t) ' 

are given by the formulas of Lemma 2.4 with Q = l^(e), and v = v{e) being an 
eigenvector of Bi{e) with the eigenvalue x\e~^ — t^^\ Similarly, [JM, Appendix A] 
shows that 

7^l 



n{o = {c-x^)Ei + no, TZ-\z) = i-Ei + 



are given by the same formulas with Q = and v being an eigenvector of Bi with 
the eigenvalue tj^\ (Note that only the off-diagonal elements of Q participate in 
the formulas.) 

Lemma 5.5. Under the assumption (5.7), we have 

AAoi?o(e)Ao ^ ^ Teo, e Aoi?i(e)Ao ^A'^ ^ Tei, e^O 
where A = e^K 

Proof. First we note that (5.7) implies that in the projective space the vector Aow(e) 
tends to 11 as e ^ 0. Next, it is easy to see that the difference Schlesinger equations 
preserve the asymptotics (5.7): 

AoBi{e;ki,...,kn)Ao^+ki-xie-^+Bi^O, e 0, (5.9) 

for any fci, . . . , A;„ e Z. In particular, 

AoCi{€)Ao^-xie-^+Bi^O, e ^ 0. 

Thus, from 

A{z, e) = Ao{e){z - Ci(e)) ■ ■ ■ {z - C„(e)) = Ao{e)z" + A,{e)z'^-' + ... 
we conclude that 

n / ^ \ 

Ai(e) = -^Cj(e) = diagonal part + A-1 ^S; + o(l) Aq. 
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By (1.3) we also know that e{sj — Si)(li(e))j,; = (^i(e))j,; lor all k ^ I. Since 
3^1 = X^^i Bi, the statement follows from the explicit formulas of Lemma 2.4. □ 

A direct computation shows that (here wc use the fact that TZqBiJZi = t^pTZoTZi 
i 1 . which follows from the explicit formulas of Lemma 2.4) 

gl = 7^o(goo + V |7^l + £;»7^l, (5.10) 

Bi = {{xi - xi)Ei + 110)31 (l~Ei+ y Z = 2,...,n. (5.11) 

V Xi- XiJ 

Let us prove (5.8) for / > 2 first. Consider the composition of the elementary 

(0 

transformation (for the difference equation) in question with <S(0, . . . , 0, 1 , 0, . . . , 0), 
see §3 for the notation. By the uniqueness part of Theorem 2.1, this is equivalent 

to making 5(0, . . . , 0, 1 , 0, . . . , 0) first, and applying the elementary transformation 
of after that. Denote the multiplier of this second elementary transformation by 
R{z, e). We have 

R{z, e) = {z- xie-^ + tf)Ei + Ao(e), R~^{z, e)=I-Ei + ^'^^^ . 

z - Xie-i +P- ' 

Using Proposition 3.6, we obtain 

{z - Biie))Riz,e) = R{z,e){z - Bi{e)). 

Substituting z = xie~^ and conjugating by Aq, we get 

AoBi{e)Ao^ - xie-^ 

= (Aoi?(a;,e-i,e)Ao^) {AoBi{e)Ao^ - xie-') {AoR-\xie-\e)Ao^) . 

Because of (5.9), the limit relations of Lemma 5.5 also hold for Rq, Ri. Using them, 
Lemma 5.5 itself, (5.7) and (5.11), we arrive at (5.8) for I > 2. 
Thus, it remains to prove (5.8) for I = 1. We have 

A{z, e) = R{z + 1, e)Ao{z - Ci(e)) ■ ■ ■ {z - Cn{e))R-\z, e) 
= Ao{z-C^{e))---{z-Cn{e)). 

The relation (5.8) for Z > 2 implies that 

AAoQ(e)Ao'A-i-a;;e-i+S,^0, e ^ 0, l = 2,...,n. 

Substituting these estimates and similar ones for C; and setting z = {w+x\)e~^ — 
t^P , we can rewrite (5.12) as follows (note that Aq is diagonal and hence it commutes 
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with Aq): 

Bn + 0{1) 



W — Xr, 



{{w -xi + e)Ei + AAoi?o(e)Ao + eB^) (^I + e ■ ■ ■ (^I + ' 

Jl-E. + fWlA^A ^ ^ (, ^ ^ .^.e- - AAoC. (e)A- A- ' 



W — Xl J V W — X\ 

\ W-X2 \ W-Xn I 



Comparing the residues of both sides &t w = x\ and looking at terms of order e, 
we see that 



e-i-AAoCi(e)Ao-iA-i^i?i, 



where Bi is given by (5.10). (Wo need to use Lemma 5.5 and the relation Ro{e)Ri{e) = 
here.) Since the difference Schlesinger equations preserve the asymptotics (5.8) 
(cf. (5.9)), we get (5.8) for I = 1, and thus for all I. 

The proof that (5.7) implies (5.8) in the case A^^^ = — Aj""-* = — 1 is very similar. 
Let us outline the necessary changes. The multipliers have the form 

R{z,e) = I-Ei+ ^iM^ ^, R-\z,e) = {z-l-Xie-'+tl'^)E,+R'o{e), 

z-1- xie-^ + t) ' 

where Ro{e), i?i(e) are constructed by Lemma 2.4 with v = v{e) being a solution 
of A*{xie-'^ - tf\e)v{e) = and Q = -Yf{e), see the proof of Theorem 2.1. 
Similarly, 

n{0 = i-Ei + -^, n-\o = {c-xi)Ei + Rl 

{ - Xl 

where TZq and TZi are as in Lemma 2.4 with v being an eigenvector of B\ with the 
eigenvalue tj^\ and Q = —yl- Similarly to Lemma 5.5, (5.7) implies 

Aoii* (e)Ao-iA-i ^ 7e* , eAAoRl{e)Ao^ ^ n\, e^O 
with A = e~^\ Similarly to (5.10), (5.11), we have 

Bi - n[ ( Soo + V ) K - KEi, 

Bi=(l-Ei + ^!^\Bi{{xi-xi)Ei+nl), l = 2,...,n. 

V Xl - Xl J 

Using the same argument with composing our elementary transformation with 

(0 

<S(0, . . . , 0, 1 , 0, . . . , 0), we prove (5.8) for I > 2. Then substituting estimates for 
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C/'s and C/'s into A{z,e) = R{z + 1, e)A{z, e)R ^{z,e), we get (with z = {w + 



W — Xi / \ W — Xi J \ w — x„ 

X ((w; - XI - e)^;, + Aoi?^(e)A^i A-i) = (/ + eBoo) 
^, (j I ^ xi6-^-AAoCi(6)A-^A-A / ^ ^ g2 + o(l) \ , + 

Comparing the residues of both sides 'At w = xi and taking terms of order e, we 
recover the estimate of type (5.8) for Ci(e), and hence for i3i(e). The proof of 
Theorem 5.3 is complete. □ 
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